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The property of isoscaling in nuclear fragmentation is studied using a simple bond percolation 
model with "isospin" added as an extra degree of freedom. It is shown analytically, first, that 
isoscaling is expected to exist in such a simple model with the only assumption of fair sampling 
with homogeneous probabilities. Second, numerical percolations of hundreds of thousands of grids 
of different sizes and with different N to Z ratios confirm this prediction with remarkable agree- 
ment. It is thus concluded that isoscaling emerges from the simple assumption of fair sampling 
with homogeneous probabilities, a requirement which, if put in the nomenclature of the minimum 
information theory, translates simply into the existence of equiprobable configurations in maximum 
, entropy states. 
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(N| 1 I- INTRODUCTION 

The experimental determination of isoscaling Q, 0, Q, § has prompted a vigorous study of its origins and its 
implications on the equation of state of asymmetric nuclear matter. Isoscaling is the property that fragment yields 

of similar, but isotopically different, reactions depend exponentially on the neutron (A) and proton (Z) numbers 

', through R 2 i{N, Z) = Y 2 (N, Z)/Yi(N, Z) « Exp[aN + f3Z], where a and /3 are fitting parameters. 

In the past, this power law expression for i?2i has been linked, under diverse approximations, to primary yields 
produced by disassembling infinite equilibrated systems in microcanonical and grand canonical ensembles [j, |6( , as 
well as in canonical ensembles 0, and it has also been observed in the framework of the grand-canonical limit of 
the statistical multifragmentation model @, in the expanding-emitting source model Q, and in the antisymmetrized 
molecular dynamics model |9j. Furthermore, under these approximations, the isoscaling parameters a and f3 have been 
found to be related to the symmetry term of the nuclear binding energy [a, , to the level of isospin equilibration 0] , 
I \ and to the values of transport coefficients pd|. 

More recently, however, it has been determined through the use of molecular dynamics Il2| that isoscaling can exist 
in purely classical systems, and that it can be created in systems fully out of equilibrium. It was also found, among 
other things, that R21 can maintain the power-law behavior even when it contains yield contributions generated at 
different times and corresponding to diverse thermodynamic conditions. 

The implications of these findings are many and very important. Isoscaling is not a quantum process; /3-decay, Pauli's 
exclusion principle and its implications for isospin selection cannot be possible causes of isoscaling. Isoscaling exists 
in finite systems out of equilibrium; expanding systems with rapidly varying temperatures and chemical potentials 
obey isoscaling. The isoscaling ratio, R21, contains contributions from different times of the reaction; its final value 
does not necessarily correspond to the thermodynamic conditions of a period of the reaction. This, of course, if not 
invalidates the scenarios and conclusions presented by previous studies, at least casts a shadow of doubt on them and, 
especially, on the assumed physical meanings of the isoscaling parameters a and (3. 

In view of the present situation, the question that needs to be addressed remains the same as in our previous 
study |12|. what produces isoscaling? Answering this question is now, in a way, simpler than in our previous work as 
many reaction variables have been eliminated out of the search. After removing the need for quantum effects and for 
thermodynamic equilibrium, it is clear that isoscaling should, then, exist in systems with little more than protons and 
neutrons without specific interactions or dynamics. The origin of isoscaling must then lie in the sampling (i.e. mode 
of fragmentation) of a conglomerate of protons and neutrons and perhaps, since we are dealing with finite systems, 
on its geometry and homogeneity. 

This work aims at elucidating the origin of isoscaling by searching for this effect in a system with the bare minimum 
number of ingredients, namely bond percolation model. After presenting the model in the next section, results of the 
percolation of hundreds of thousands of three-dimensional grids are presented in section IIIII followed with a number 
of conclusions in section llVl 
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II. THE PERCOLATION MODEL 



In order to explore the behavior of isoscaling we use a three dimensional bond percolation model. This model was 
first applied in nuclear multifragmentation by Bauer et al. ^| an d used by many groups 0, 0, 0, 0] ever since. 
In the usual bond percolation model, a fragmenting nucleus is represented by a three-dimensional cubic lattice, and 
individual "nucleons" by nodes on the lattice. All nodes start with bonds to all nearest-neighbors, which represent 
nucleon-nucleon interactions. These bonds are then attempted to be broken statistically according to a probability 6, 
thus producing clusters of connected nodes which are interpreted as fragments. 

To this usual model, we add the "isospin" of the nodes as an extra degree of freedom. With this method, lattices 
with different ratios of "protons" and "neutrons" can be constructed and ruptured, producing cluster yields which 
can then be used to construct the ratio i?2i- Besides the usual assumptions of bond percolation, the following three 
conditions are added in this analysis: i) the bond breaking probability, 6, is spatially homogeneous and identical for 
pp, pn, and nn bonds, ii) the probability of a node for having isospin "up" (p) or "down" (n) is spatially homogeneous, 
and Hi) the number of protons and neutrons is fixed from the beginning. 



A. Geometrical arguments leading to isoscaling 



In this model, fragments are obtained when bonds are broken with a given probability b. If we consider the infinite 
size limit is customary to express the number of fragments as the number of fragments per node 
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where L 3 is the size of the lattice measured in nodes, N A is the number of fragments of size A, and tia the number 
of fragments of size A per node. This last quantity can be written in the following way 



n A 



gAatil-byb 1 



(2) 



where t stands for the perimeter of the cluster (number of bonds to be broken in order to isolate the cluster composed 
by A nodes), and a is the number of bonds linking the A nodes, g Aa t is the number of cluster configurations with 
size A, perimeter t and a activated bonds. As an illustration, figure shows the different terms appearing in the 
expression for 714 in the two dimensional case. The resulting expression for this term is 



n 4 = 14(i - bfb 10 + 4(1 - b) 3 b 9 + (I - b) 4 b 8 



(3) 



To include the isospin degree of freedom, a given node belonging to a cluster of size A will be considered to be a 
proton with probability p. Then 
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with az being the number of ways of building a cluster with Z protons and N = A — Z neutrons: az = Al/NlZl The 
number of fragments per node with A particles and Z protons is then defined as: 



n>A,z = n A [azp Z q^ A Z) ] = n A a z p Z q N ■ 
Focusing now on the isoscaling problem, the quantity of interest is the quotient: 

Y 2 (A,Z) 



R21 = 



Yi(A,Z) 



(5) 
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With Y2 representing the yield for the reaction involving the neutron rich nuclei. 
For the percolation model, R21 takes the form 



[nA,zh = P2I2 
[nA,z]i pfgf \pi 



R21 = 
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FIG. 1: Cluster structures in two dimensional percolation. In this figure we show the 6 possible realizations of clusters of size 
4 in a two dimensional grid. Circles denote nodes, straight lines denote active links. Together with each diagram we show the 
corresponding term in the sum m. 



in which all the geometrical and combinatorial terms cancel out and only the part related to the occupancy probabilities 
remains. In this way 



i?2i = exp < In 



N' 
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Nhi | — 
<7i 



Pi 



exp (aN + (3Z) 



(8) 



Which is the standard expression of the isoscaling coefficient, but now with a clear microscopic interpretation for the 
isoscaling parameters a and [3. 

For finite lattices, it is convenient to study R21 calculated in terms of the number of fragments Na,z instead of 
the above derivation in terms of the number of fragments per node, ua,z- This can be achieved by using the above 
derived expressions for the constants a and f3. Equation Q can be rewritten as 
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where the probabilities have been approximated by pi ~ Z/A{. This immediately yields the usual expression: 

R21 = Cexp(a7V + /3Z) (10) 

with a = ln(q 2 /qi), (3 = ln(p 2 /Pi) and C = exp[ln(pi/p 2 )] = P1/P2 = A 2 /A x 

It should be kept in mind that the main assumptions in this derivation are that the probability p is homogenous 
and that the number of fragments per node (per particle) can be approximated by the corresponding infinite size 
limit. We now turn to a numerical verification of these predictions. 

III. NUMERICAL EXPERIMENTS 

To explore the isoscaling phenomena in the framework of the above defined percolation model, two different types 
of calculations were performed: fixing the number of protons in the lattice, and using a fixed probability to assign the 
isospin to the nodes. 

A. Fixed proton number 

In this case, the R 2 \ was obtained using the yields of two lattices, one of size 6x6x6, i.e. with A = 216 nodes, 
and a second one of 7 x 7 x 7 with A = 343 nodes, both with the number of protons fixed to 108. These grids, had, 
then, Z = N = 108, and probabilities p x = 0.5 and q 1 = (1 - p x ) = 0.5 for the 6 x 6 x 6 grid, and Z = 108, N = 235, 
p 2 = 0.315 and q 2 = (1 - p 2 ) = 0.685, for the 7 x 7 x 7 grid. 

For these probabilities, the isoscaling coefficients are a = ln(0. 6852/0. 5) = 0.315 and (3 — ln(0. 3148/0. 5) = —0.463, 
and the constant C = pi/p 2 = 0.5/0.3148 = 1. 588 3. The results of 200, 000 percolations of each of these grids are 
displayed in panel a) of figure 

B. Fixed p and n occupancy probability 

In this second case, the two isotopically different grids were constructed with the same sizes for both, but with two 
different protons and neutrons occupation probabilities. Several cases were studied. 

First, 6x6x6 grids were constructed using p\ = 0.5 and q\ — (1 —px) = 0.5 which yields a total number of protons 
(in average) of (Z) = 108 and (N) — 108, and p 2 = 0.33 and q 2 = (1 — p 2 ) = 0.67 which yields a total number of 
protons (in average) of (Z) = 72 and (N) = 144. In this case the coefficients are : a = ln(0.67/0.5) = 0.292 67, and 
(3 = ln(0.33/0.5) = —0.415 52. The results obtained from this numerical exercise are displayed in panel b) of figure|2J 
the fact that isoscaling is well reproduced by this purely geometrical model is obvious from this figure. 

Next, to investigate size effects, a similar case was constructed using a smaller lattice of size 5x5x5 with the 
same probabilities as in the previous case. The corresponding results are displayed in panel c) of figure [21 again the 
property of isoscaling is apparent. 

Finally, to investigate the effect of different occupation probabilities on R 2 \, one more case was studied. In this 
case, a lattice of 6 x 6 x 6 was populated with protons and neutrons according to the probabilities: p\ — 0.5 and 
q\ = (1 — px) = 0.5. For the neutron rich partner the probabilities used were p 2 = 0.42 and q 2 = (1 — p 2 ) = 0.58 
which yields a total number of protons (on average) of (Z) = 91 and (N) = 125. In this case the coefficients are 
a = ln(0. 58/0.5) = 0.148 42, and (3 = ln(0.42/0.5) = -0.174 35. The corresponding results are shown in panel d) of 
figure |3 

For each of these cases 200, 000 configurations were generated. The bond breaking probability was chosen as b = 0.31 
in order to get good statistics in an ample range of masses. [Attempts with different values of the breaking probability 
b demonstrated that the ratio R 2 \ is independent of b.] The figure also shows a comparison of the numerical results 
to the theoretical predictions of section II I II lines denote theoretical predictions whereas symbols denote numerical 
simulations. The agreement between the theoretical predictions and the results of the simulations is remarkable. 

IV. CONCLUSIONS 

We have studied the isoscaling phenomenon in the frame of percolation model. We have derived exact analytic 
expressions for the infinite case and approximate ones for the finite case. We have performed numerical simulations for 
not too big systems (125, 216 and 343 "particles") with different relative populations of N:Z. The excellent agreement 
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FIG. 2: Isoscaling in percolation. Comparison of theoretical predictions for isoscaling in the nuclear percolation model and the 
corresponding results from numerical simulations. Lines denote theoretical predictions, symbols denote numerical simulations 



between numerical simulations and theory indicate that isoscaling emerges from the simple assumption of fair sampling 
with homogeneous probabilities. On the other hand, this property can be seen as a minimum information approach, 
i.e. all configurations are equiprobable, as such this analysis can be interpreted in the frame of a maximum entropy 
approach. This indicates that the information about effects due, for example, to the asymmetry term in the equation 
of state, is in the absolute values of the parameters a and j3, and not in the isoscaling property itself. 
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